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Abstract 

We present existence, uniqueness and continuous dependence results 
for some kinetic equations motivated by models for the collective behavior 
of large groups of individuals. Models of this kind have been recently pro- 
posed to study the behavior of large groups of animals, such as flocks of 
birds, swarms, or schools of fish. Our aim is to give a well-posedness theory 
for general models which possibly include a variety of effects: an interac- 
tion through a potential, such as a short-range repulsion and long-range 
attraction; a velocity-averaging effect where individuals try to adapt their 
own velocity to that of other individuals in their surroundings; and self- 
propulsion effects, which take into account effects on one individual that 
are independent of the others. We develop our theory in a space of mea- 
sures, using mass transportation distances. As consequences of our theory 
we show also the convergence of particle systems to their corresponding 
kinetic equations, and the local-in-time convergence to the hydrodynamic 
limit for one of the models. 
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1 Introduction 



The description of the collective motion (swarming) of multi-agent aggregates 
resulting into large-scale structures is a striking phenomena, as illustrated by the 
examples provided by birds, fish, bees or ants. Explaining the emergence of these 
coordinated movements in terms of microscopic decisions of each individual 
member of a swarm is a hot matter of research in the natural sciences [H [TTJ EH] . 
The formation of swarms and milling or flocking patterns have been reported 
in animals with highly developed social organization like insects (locusts, bees, 
ants, ...) [11], fishes [UE] and birds (6] EH] but also in micro-organisms as myxo- 
bacteria [22]. Moreover, the understanding of natural swarms has been used as 
an engineering design principle for unmanned artificial robots operation [5j [30] . 

The physics and applied mathematics literature has proliferated and sprung 
in this direction in the recent years trying to model these phenomena, mainly 
based on two strategies of description: individual-based models or particle dy- 
namics [351 Ell [23J El [261 [HI [III [HI H2J H3l [Ml ES] and continuum models 
based on PDEs for the density or for the momentum of the particle ensem- 
ble EH1 [Ml 1311 EH HI- The key feature to explain is the emergence of self- 
organization: flocking, milling, double milling patterns or other coherent be- 
havior. 

Particle descriptions usually include three basic mechanisms in different re- 
gions: short-range repulsion zone, long-range attraction zone and alignment 
or orientation zone, leading to the so-called three-zone models. In addition, 
some of them incorporate a mechanism for establishing a fixed asymptotic 
speed/ velocity vector of agents, as is usually observed in nature. Some of the 
models only consider the orientation vector and not the speed in their discrete 
version. The main differences of all these models reside in how these three in- 
teractions are specifically considered. We will mainly work with two generic 
examples in which several of the effects above are included, namely the model 
for self-propelled interacting particles introduced by D'Orsogna et al in [15] and 
the model of alignment proposed by Cucker and Smale [321 [13] . 

Together with particle and continuum models based on macroscopic den- 
sities, there has been a very recent trend of mesoscopic models by means of 
kinetic equations for swarming EH1 [71 [TH1 [5] . In these models one works with 
a statistical description of the interacting agent system. Let us represent by 
x £ R d the position, where d > 1 stands for the physical space dimension, 
and by v £ M. d the velocity. We are interested in studying the evolution of 
/ = f(t, x, v) representing the probability measure/density of individuals at po- 
sition x, with velocity v, and at time t > 0. These are the models we study in 
the present paper. Given that we cover a variety of them, we refer the reader 
to Section [2] for a more detailed presentation of the equations. 

These kinetic models bridge the particle description of swarming to the hy- 
drodynamic one as already discussed in [20j [Jj [19] . The main key idea is that 
solutions to particle systems are in fact atomic-measure solutions for the ki- 
netic equations, and solutions to the hydrodynamic equations are solutions of a 
special form to the kinetic equation; see Section [5.21 for more details. 
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In some cases, suitable compactness arguments based on the stability prop- 
erties in distances between probability measures allow to construct a well- 
posedness theory for a kinetic equation. Such an approach was done for the 
Vlasov equation in classical kinetic theory [37J 111 Ell EH] with several nice reviews 
in [551 1321 E] ■ Of these references, [T3] uses the Monge-Kantorovich- Rubinstein 
distance (the one we use in the present paper); the others, as well as the recent 
work |19j for the kinetic Cucker-Smale model, use an approach based on the 
bounded Lipschitz distance. 

In this paper we present a generic approach to the well-posedness of many 
of these models in the set of probability measures in phase space based on the 
modern theory of optimal transport [36 . In fact, we will use the well-known 
Monge-Kantorovich-Rubinstcin distance between probability measures instead 
of the bounded Lipschitz distance. Its better duality properties actually make 
this technical approach easier in terms of estimates leading to one of our crucial 
results: a stability property of solutions to swarming equations under quite 
general conditions. 

We derive some consequences from this stability estimate. First, we prove 
the mean-field limit, or convergence of the particle method towards a measure 
solution of the kinetic equation. This mean field limit is then established without 
any resorting to the BBGKY hierarchy or the molecular chaos hypothesis [3, 
[TJQIJ]. Second, we show the stability for arbitrary times of the hydrodynamic 
solutions, assuming they exist, although with constants depending on time. 
Finally, the stability result can be used to obtain qualitative properties of the 
measure solutions of the kinetic equations, as it has been done in [8J for the 
kinetic Cucker-Smale model. 

This strategy is quite general, and we first demonstrate its use in a particular 
kinetic model introduced in [7] for dealing with the mesoscopic description and 
certain patterns not covered by the particle model proposed in |15j . Other 
models are treated by the same procedure in subsequent sections, as the kinetic 
Cucker-Smale model proposed in [20] for the original alignment mechanism in 
[T21 [TB"] . the models studied in [2H US], or any linear combination of these 
mechanisms. We finally give general conditions on a model that are sufficient 
for our well-posedness results to be valid. 

Let us comment on some limitations of the method we use. The first one is 
that, as we work with solutions in a weak measure sense, we have to require our 
interaction terms to be locally Lipschitz in order in order to carry out the theory. 
This is a well-known limitation in the literature for working with the mean field 
limit and measure solutions, see [3H [3T] and the references therein. A less 
fundamental one is that we always work with compactly supported solutions. 
One could probably develop a theory substituting this condition by a suitable 
control on moments of the solution, and then adapting the estimates to this 
setting; however, in the present paper we do not pursue further extensions in 
this direction. 

Next section does a simple and brief review of the main interacting particle 
systems under analysis and the needed concepts for the Monge-Kantorovich- 
Rubinstein distance between probability measures. The third section is devoted 
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to the proof of the main result of existence, uniqueness and stability of measure 
solutions to the particular swarming equations introduced in [7]. Section 4 
generalizes this approach to a general family of these equations. Finally, section 
5 draws some consequences of the stability property: the convergence of the 
particle method and the mean-field limit are proved for the general model, while 
the stability in a finite time interval of hydrodynamic solutions is shown for the 
swarming model used in Section [3l 



2 Preliminaries 

In this section we introduce the models mentioned in the introduction. We give 
some particular representative cases and specify the models to which our results 
apply. Also, we recall some notions about optimal transport that shall come in 
handy. 

2.1 Main Kinetic Models 

The particle model proposed in [15] reads as: 

^f = vu (i = l,...,N) 

^ = (a-p \ Vi \ 2 ) Vi - 1 VU(\ Xi ~ Xj \), (i = l,...,N). 

where a, ft are nonnegative parameters, U : R d — > K is a given potential 
modeling the short-range repulsion and long-range attraction typical in these 
models, and N is the number of particles. Here, the potential has been scaled 
depending on the mass of each particle as in [7] , where we refer for further dis- 
cussion. The term corresponding to a models the self-propulsion of individuals, 
whereas the term corresponding to /3 is the friction assumed to follow Rayleigh's 
law. The balance of these two terms imposes an asymptotic speed to the agent 
(if other effects are ignored), but does not influence the orientation vector. A 
typical choice for U is the Morse potential which is radial and given by 

U(x) = k(\x\) with fc(r) = -C A e~ r,tA + C R e~ r/eR , 

where Ca, C r and £a,£r are the strengths and the typical lengths of attraction 
and repulsion, respectively. This potential does not satisfy the smoothness 
assumption in our main theorems but the qualitative behavior of the particle 
system does not depend on this particular fact [15j . In fact, a typical potential 
satisfying all of our hypotheses is 

The kinetic equation associated to this particle model as discussed in [7] 
gives the evolution of / = fit, x, v) as 

dtf + v- V x f - iVU * p) ■ W v f + div„((a - p \v\ 2 )vf) = 0, (1) 
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where p represents the macroscopic density of /: 



p(t,x) := / f(t,x,v)dv for (>0,i£l li . (2) 



In the Cucker-Smale model, introduced in [HI [13], the only mechanism taken 
into account is the reorientation interaction between agents. Each agent in the 
swarm tries to mimick other individuals by adjusting/ averaging their relative 
velocity with all the others. This averaging is weighted in such a way that 
closer individuals have more influence than further ones. For a system with TV 
individuals the Cucker-Smale model reads as 



N 



dxi 




~dt 




dvt 


1 


dt 


~ N 



with the communication rate w(x) given by: 



-3\,- ( 1 + | a .._ a . j |2)7 

for some 7 > 0. This particle model leads to the following kinetic model [20l 

nan]: 

|+«-VJ = V»-[([/]/] (3) 

where £[f](x,v,t) = (H * /) (x,v,t), with H(x,v) — w(x)v and * standing for 
the convolution in both position and velocity (x and v). We refer to [12l H3l 18] 
for further discussion about this model and qualitative properties. 

Moreover, quite general models incorporating the three effects previously 
discussed have been considered in [231 ESI ■ I R particular, they consider that N 
individuals follow the system: 



(4) 

v az 

where F^ is the self-propulsion autonomously generated by the ith-individual, 
while F? is due to interaction with the others. The model in Section 3 corre- 
sponds to F^ = (a ~ (3 \vi\ 2 )vi, while the term Ff- — —(3vi is considered in [2"3"] . 
and Ff- = a, — j3vi in [3U[25]. Here, dj is an autonomous self-propulsion force 
generated by the ith-particle, and may depend on environmental influences and 
the location of the particle in the school. The interaction with other individuals 
can be generally modeled as: 

N N 




n - f[- + f^ = y: 9± {\xi - ^d^tti + e mi* 

J= l 3=1 
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Here, g+ and h + (<?_ and hJ) are chosen when the influence comes from the 
front (behind), i.e., if (xj — Xi) • Vi > (< 0); choosing g + ^ g_ and h + ^ h- 
means that the forces from particles in front and those from particles behind 
are different. The sign of the functions g±(r) encodes the short-range repulsion 
and long-range attraction for particles in front of (+) and behind (-) the ith- 
particle. Similarly, h + > (< 0) implies that the velocity-dependent force 
makes the velocity of particle i get closer to (away from) that of particle j. 

In the next sections we will be concerned with the well-posedness for measure 
solutions to JT]) , ([3J and generalized kinetic equations including the correspond- 
ing to the TV- individuals model in (QJ. 



2.2 Preliminaries on mass transportation and notation 

Let us recall some notation and known results about mass transportation that 
we will use in the next sections. For a more detailed approach, the interested 
reader can refer to [9j [36] . 

We consider the space of probability measures Vi (R d ) , consisting of all prob- 
ability measures on R d with finite first moment. In Vi(Mr) a natural concept of 
distance to work with is the so-called Monge-Kantorovich-Rubinstein distance, 

Wx(f,g) = sup | J <p(P)(f(P)-g(P))dP , <p € Lip(R d ), Lip(^) < lj , (5) 

where Lip(R d ) denotes the set of Lipschitz functions on R d and Lvp(<p) the 
Lipschitz constant of a function if. Denoting by A the set of transference plans 
between the measures / and g, i.e., probability measures in the product space 
R d x M. d with first and second marginals / and g respectively, then we have 

Wx{f,g) = inf ( / |P X - P a | d^{P u P 2 )\ (6) 

by Kantorovich duality. 'Pi(K d ) endowed with this distance is a complete metric 
space. In the following proposition we recall some of its properties. We refer to 
[55] for a survey of these basic facts. 

Proposition 2.1 (Wi -properties). The following properties of the distance W\ 
hold: 

i) Optimal transference plan: The infimum in the definition of the dis- 
tance W\ is achieved. Any joint probability measure II satisfying: 

Wx(f,g)= f \P 1 -P 2 \dn o (P 1 ,P 2 ). 

JS. d xM d 

is called an optimal transference plan and it is generically non unique for 
the W\-distance. 

ii) Convergence of measures: Given {fk}k>i and f in Vi(M. ), the fol- 
lowing three assertions are equivalent: 



6 



a) W\(fk,f) tends to as n goes to infinity. 

b) fk tends to f weakly-* as measures as k goes to infinity and 



sup / \v\ fk{v) dv — » as i? — > +oo. 

fc>l J\v\>R 



fc>l J\v\>R 

c) fk tends to f weakly-* as measures and 

v\fk(v)dv^ / \v\f(y)dv as n^+oo. 

Throughout the paper we will denote the integral of a function tp = ip(x) 
with respect to a measure p by J p{x)p{x) dx, even if the measure is not abso- 
lutely continuous with respect to Lebesgue measure, and hence does not have 
an associated density. 

Given a probability measure / 6 V\ (R d x R d ) we always denote by p its first 
marginal, written as follows by an abuse of notation: 



P( x ) '■= / f{x,v)dv. (7) 
To be more precise, p is given by its action on a C° function <f> : M. d — s- M, 
p(x)(f>(x) dx = / f(x,v)(p(x)dxdv. 



For T > and a function / : [0,T] -> x M d ), it is understood that p is 

the function p : [0,T] — > ^(R ) obtained by taking the first marginal at each 
time t. Whenever we need to indicate explicitly the dependence of p on /, we 
write p[f] instead of just p. 

We denote by Br the closed ball with center and radius R > in the 
Euclidean space R™ of some dimension n. When we need to explicitly indicate 
the dimension of the space, we will write B R . For a function H : R" — >• ]R m , we 
will write LipR(H) to denote the Lipschitz constant of H in the ball Br C R". 
For T > and a function H : [0,T] x R" -> R m , H = H(t,x), we again write 
LipR(H) to denote the Lipschitz constant with respect to x of H in the ball 
-Br C R n ; this is, LipR(H) is the smallest constant such that 

|fl'(t,a;i)-fl'(t,xa)| < Lip R (H) \ Xl - x 2 \ for aU xi,x 2 e B R ,t e [0,1*]. 

For any such function H , we will denote the function depending on i at a fixed 
time t by this is, H t (x) := ff(t, x). 



3 Well-posedness for a system with interaction 
and self-propulsion 

In this section we consider eq. ([1]). In this model (and in fact, in every model 
considered in this paper) the total mass is preserved, and by rescaling the equa- 
tion and adapting the parameters suitably one easily sees that we can normalize 
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the equation and consider only solutions with total mass 1. We will do so and 
reduce ourselves to work with probability measures. 

3.1 Notion of solution 

In order to motivate our definition of solution to equation ([1} let us consider 
for a moment a general field E instead of —VU * p. Precisely, fix T > and a 
function E : [0,T] xR d -^M. d such that: 

Hypothesis 3.1 (Conditions on E). 1. E is continuous on [0,T] x M. d , 

2. For some Ce > 0, 

\E(t,x)\<C E (l + \x\), for allt,x e[0,T]xR d , and (8) 

3. E is locally Lipschitz with respect to x, i.e., for any compact set K C M. d 
there is some > such that 

\E(t,x)-E(t,y)\<L K \x-y\, t € [0,T\, x,y G K. (9) 

We consider the equation 

dtf + v ■ Vje/ + E ■ V„/ + div w ((a - [3 \v\ 2 )vf) = 0, (10) 

which is a linear first-order equation. The associated characteristic system of 
ode's is 

V, (11a) 

E{t,X) + V(a-/3\V\ 2 ). (lib) 

Lemma 3.2 (Flow Map). Take afield E : [0, T] xM d -)• R d satisfying Hypothesis 
\3.1\ Given (Xo, Vq) € M. d x M. d there exists a unique solution (X, V) to equations 
ifria P -lfTlb )) in C 1 ([0,T];E d x R d ) satisfying X(0) = X and V(0) = Vq. In 
addition, there exists a constant C which depends only on T, \Xq\, \Vq\, a, fj 
and the constant Ce in eq. ((HJ, such that 

\{X{t),V(t))\ < \(X Q ,V )\e ct for alltG [0,T]. (12) 

Proof. As the field E satisfies the regularity and growth conditions in Hypothesis 
13 - 1 1 standard results in ordinary differential equations show that for each initial 
condition (X(0), V(0)) G K d x K d this system has a unique solution defined on 
[0, T) (the only term in the equations which does not grow linearly is — j3V \V\ 2 , 
and it makes \V\ decrease, so the solution is globally defined in time). The 
bound (fT2|) on the solutions follows from direct estimates on the equation, using 
the linear growth of the field E. □ 



dt 
d 

dt 



-X 



V 
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Calling P = (X, V), the system (|llgi[) - (|l lbt) can be conveniently written as 

j t P = * E (t,P), (13) 

where ^ E ■ [0, T] x R d x R d -> R d x M d is the right hand side of eqs. (fTTa|) , 
(jllbp . When the field -E is understood we will just write *f? instead of ^ E . Using 
this notation, equation (|10p can also be rewriten as 

df 

-± + div(* B /) = 0. (14) 
We can thus consider the flow at time t £ [0, T) of eqs. (fTTj) . 
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d x R d -» R d x R d . 



Again by basic results in ode's, the map (t,x,v) H> T E {x,v) — (X, V) with 
(X, V) the solution at time t to (fl2|) with initial data (x, v), is jointly continuous 
in (t, w). For a measure jo e Pi(M d x R d ) it is well-known that the function 

f :[0,T)^P 1 (R d xR d ), t ->/*:= 72#/„ 



is a measure solution to eq. p0|) . i.e., a solution in the distributional sense. Here 
we are using the mass transportation notation of push-forward: f t — T^#/o is 
defined by 

((x,v) f(t,x,v)d(x,v) = ((T E (x,v)) f Q (x,v)d(x,v), (15) 

for all C G C£(R 2d ). Note that in the case where the initial condition /o is regular 
(say, C£°) this is just a way to rewrite the solution of the equation through the 
method of characteristics. This motivates the following definition: 

Definition 3.3 (Notion of Solution). Take a potential U € C 1 (M ti ) such that 

|Vl/(ar)| < C7(l + |ar|), x G R d , (16) 

for some constant C > 0. Take also a measure fo S Pi (II x R d ), and T E 
(0, oo]. We saj/ t/iat a function f : [0,T] — > Vi(R d x M d ) is a solution of the 
swarming equation JT]) iwii/i initial condition fo when: 



1. The field E[f] — —VU * p satisfies the conditions in Hvvothesis \3.1[ 

2. It holds f t = T% m #fo. 

Remark 3.4. This definition gives a convenient condition on U so that a mea- 
sure solution in V\{R d x R d ) makes sense. One can weaken the requirement on 
U in this definition as long as the requirements on f are suitably strengthened 
(e.g., one can allow a faster growth of the potential if one imposes a faster decay 
of f , or less local regularity of U if one assumes more regularity of f), but we 
will not consider these modifications in the present paper. 
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Since we ask the gradient of the potential to be locally Lipschitz, we cannot 
consider potentials with a singularity at the origin. This is a strong limitation of 
the classical theory, and is considered a difficult problem for the mean-field limit. 
As for the existence theory, if one wants to consider more singular potentials, 
one can work with functions f which are more regular than just measures, so 
that X7U * p becomes locally Lipschitz and a parallel existence theory can be 
developed. 

3.2 Estimates on the characteristics 

We gather in this section some estimates on solutions to the characteristic equa- 
tions CLE]). In this section we fix T > and fields E, E X ,E 2 : [0,T] x R d -> W l 
which are assumed to satisfy Hypothesis 13.11 and we consider their correspond- 
ing characteristic equations (ITTj) . Recall that f e is a shorthand for the right 
hand side of (1TT1) . as in (|T3|) . 

We first gather some basic regularity results for the function which defines 
the right hand side of eqs. (|lla[) - (|llb|) : 

Lemma 3.5 (Regularity of the characteristic equations). Take a field E : 
[0, T] x W 1 — )> M. d which satisfies Hypothesis lff.il Consider a number R > 
and the closed ball B R C M d x R d . 

1. $ e is bounded in compact sets: For P = (X, V) £ Br and t £ [0, T], 

\*E(t,P)\<C 

for some C > which depends only on a, f3, R, and \\E\\ Ltx> ,r T y xBlj y 

2. $ e is locally Lipschitz with respect to x,v: For all P\ = (Xi,V\), P2 = 
(X 2 ,V 2 ) m B R , andt£ [0,T], 

|tf B (t, Pi) - tt B (t, P 2 )| < C(l + Li PR {E t )) \Pi-P 2 \, 

for some number C > which depends only on a and fj. 

Proof. This can be obtained by a direct calculation from eqs. (II la[) (|1 11d|) . □ 

Lemma 3.6 (Dependence of the characteristic equations on E). Take two fields 
E 1 , E 2 : [0,T] x R d — > M. d satisfying Hupothesis \3 . 1\ and consider the functions 
^ e 1 , ^ e 2 which define the characteristic equations (|11|) as in eq. (|13j) . Then, 
for any compact (in fact, any measurable) set B, 

ll^-v^H^ < \\E l - E 2 \\ LOO{By 

Proof. Trivial from the expression of ^e x i ^e 2 - d 

Now we explicitly state some results which give a quantitative bound on the 
regularity of the flow T E , and its dependence on the field E. 
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Lemma 3.7 (Dependence of characteristics on E). Take two fields E X ,E 2 : 
[0,T] x W l ->• R d satisfying Hypothesis \3J\ and a point P° € M d x R d . Take 
R > 0, and assume that 

|T|i(P°)| <R, \T^2(P°)\<R forte[0,T]. 

Then for t G [0, T] it holds that for some constant C which depends only on a, 
(3, R andUpxiE 1 ) 

p Ct _ i 

\Tei(P°) ~ T~e^(.P°)\ < sup H^-^ll (17) 

° se[o,T) ' 

Proof. For ease of notation, write P,(t) = 7^(P°) = (Xi(t),Vi(t)), for i = 1, 2, 
t G [0,T]. These functions satisfy the characteristic equations (| 1 lj) : 

j t P i = V E4 (t,P i ), P(0) = P°, fori = 1,2. 



Then, for i G [0,P], and using Lemmas 13.51 and 13.61 

\Pi{t) - P 2 {t)\ < f \^ El (s,P 1 (s)) -* E2 (s,P 2 (s))\ ds 
Jo 

< f \* E i(s,Pi(s))-* E i(s,P 2 (s))\ds 
Jo 

+ f \^ E i(s,P 2 {s))-^ E2 (s,P 2 ( S ))\ds 

Jo 

<C [ |P 1 (s)-P 2 (s)| ds+ [ \\El-E, 
Jo Jo 



2 11 

s \\L°°(B R ) 



ds 



where C is the constant in point 2 of Lemma l3.51 which depends on a, (3, R and 
the Lipschitz constant of E 1 with respect to x in the ball Br. By Gronwall's 
Lemma, 



\Px(t) - P 2 (t)\ < f e c ^\\El-E 
Jo 



2 II ds 



e ct - 1 

^—7^ SU P \\ E . -^L~ (Br) , 

which finishes the proof. □ 

Lemma 3.8 (Regularity of characteristics with respect to initial conditions). 
Take T > and a field E : [0,T] x M d ->• M d satisfying Hypothesis [OJ Tate 
a/so Pi,P 2 G K d x K d and P > 0, and assume that 

\Te(Pi)\<R, \Te(P 2 )\<R te[0,T]. 

Then it holds that 

\Te{Px) - 7l(P 2 )| < | Pi -P 2 |e c ^( Li P«(^)+i)^ ; ; G [0,T], (18) 
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for some constant C which depends only on R, a and ft. Said otherwise, is 
Lipschitz on Br C R d x K d , with constant 

Li Pfl (Ti) < e c£(u Pa (E.)+i)* t g [0)T] _ 

Proo/. Write P 4 (f) = T|(Pi) = P$(t)), for i = 1,2, f e [0, T]. These 

functions satisfy the characteristic equations (fTTj) : 

^Pi = * B (t,Pi), Pi(0)=Pi, for i= 1,2. 
di 

For t € [0,T], using Lemma 1331 

|Pi(t)-P 2 (t)| < |Pi-f 2 |+ / |*b(s,Pi(s))-*b(s,P 2 (s))Ms 

Jo 

< |Pi-P 2 |+C / (Up R (E s ) + l)\P 1 (s)-P 2 (s)\ ds 
Jo 

We get our result by applying Gronwall's Lemma to this inequality. □ 

Lemma 3.9 (Regularity of characteristics with respect to time). Take T > 
and a field E : [0, T] x R d -» R d satisfying Hypothesis\3Ji Take P° G R d x M d ; 
i? > and assume that 

\T*(P°)\<R, t€[0,T}. 

Then it holds that 

\r E (P Q )-T£(P°)\<C\t-s\ for s,t G [0,T], (19) 

/or some constant C which depends only on a, /3, R and \\E\\ Laa ,r Q r] X £ B y 

Proof. Direct by definition of T^(P°) and from point 1 of Lemma [3~5l as we are 
assuming that 7e(P°) remains on a certain compact subset of R d x M. d . □ 

3.3 Existence and uniqueness 

Theorem 3.10 (Existence and uniqueness of measure solutions). Take a po- 
tential U G C 1 (M d ) such that VU is locally Lipschitz and such that for some 
C>0, 

\VU(x)\ < (7(1 + forallxeR d , (20) 

and fo G Vi(M. d x M d ) with compact support. There exists a solution f on 
[0, +oo) to equation ([TJ with initial condition fo in the sense of Definition 1 3. ffl 
In addition, 

f G C([0, +oo); Pi(R d x R d )) (21) 
and t/iere is some increasing function R = R(T) such that for all T > 0, 

supp f t C B fl(T) C R d x R d /or a// i G [0, T]. (22) 

This solution is unique among the family of solutions satisfying (|21[) and (|22j) . 
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The rest of this section is dedicated to the proof of this result, for which 
we will need some previous lemmas. We begin with a general result on the 
transportation of a measure by two different functions: 

Lemma 3.11. Let P\,P 2 : W l — > R d be two Borel measurable functions. Also, 
take f eVi{R d ). Then, 

W 1 (P 1 #/,P 2 #/) < ||P X - P 2 |U (supp/) • (23) 

Proof. We consider a transference plan defined by n := (Pi x P 2 )#f. One can 
check that this measure has marginals -Pi#/, i"2#/- Then, 

Wi(Pi#/,P 2 #/)< I \x-y\Tr(x,y)dxdy 

JR d xR d 

IP^x) - P 2 (x)\ f(x)dx<\\P 1 -P 2 \\ L ^ {Bnppf) , 

which proves the lemma. □ 

Lemma 3.12 (Continuity with respect to time). Take T > and a field E : 
[0, T] x M. d — > M. d in the conditions of Hypothesis HOI Take also a measure f on 
M. d x M. d with compact support contained in the ball Br. 

Then, there exists C > depending only on a, f3, R and \\E\\ LOO ,r T -i x g R \ 
such that 

Wi(Ti#f, T%*f) <C\t-s\, for any t, s 6 [0, T]. 

Proof. From Lemma 13.111 and the continuity of characteristics with respect to 
time, Lemma T3.91 we get 



Wi(T£#f,n#f) < l|TJ-7!||^ (supp/ ) <C\t- s \, 
for some C > which depends only on the quantities in the lemma. □ 

Lemma 3.13. Take a locally Lipschitz map T : R d -)• M d and f,g £ Vi(R d ), 
both with compact support contained in the ball Br . Then, 

Wi(T#f, T#g) < L W t (f, g), (24) 

where L is the Lipschitz constant of T on the ball Br. 

Proof. Set 7r to be an optimal transportation plan between / and g. The mea- 
sure 7 = (T x T)#7r has marginals T#/ and T#g, as can be easily checked, 
so we can use it to bound W\(T#f,Tifg). 

Wi(T#f, T4^g) < / \z — w\ j(z, w) dzdw = / \T{z) — T{w)\ tt(z, w) dz dw 

jR d xR d Jm d xR d 

<L \z — w\ w(z, w) dzdw = L Wi(f, g), 

JR d xR d 

using that the support of it is contained in Br x Br, as both / and g have 
support inside Br. □ 
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Recalling that E[f] := VU*p, the properties of convolution lead immediately 
to the following information: 

Lemma 3.14. Take a potential U : R d — > K in the conditions of Theorem ] 3. 1(A 
and a measure f € 7 3 i(K d X R d ) with support contained in a ball B R . Then, 



\\E[f]\\ L ^ {BR) < ||VC/|| L „ (Bafl) , (25) 

and 

Up R (E[f}) < Up 2R (VU). (26) 
Lemma 3.15. For f,g G Vi(R d x R d ) and R > it holds that 

\\E[f] - E[g]\\L°°(B R ) < Up 2R (VU)W 1 (f,g). (27) 

Proof. Take tt to be an optimal transportation plan between the measures / 
and g. Then, for any x £ B R , using that n has marginals / and g, 



E[f](x) - E[g](x) = / (p[f](y) - p[g](y))VU(x - y) dy 

f(y,v)\7U(x — y)dydv— / g(z,w)'VU(x — z) dz duo 

JR d xS. d 

(VU(x — y) — VU(x — z)) dir(y, v, z, w). 



/R4d 

Taking absolute value, 



\E[f](z) - E[g]{z)\ < f \VU(x-y)-VU(x-z)\dir(y,v,z,w) 

< Lip 27? (VC/) f \y-z\ dir{y,v,z,w) <Lip 2R (VU)Wi(f,g), 

JWL 4d 

using that n(y, v, z, w) has support on B R x B R C R 4d . □ 

We can now give the proof of the existence and uniqueness result. 

Proof of theorem WUA Take / € T jl (R d x R d ) with support contained in a ball 
B R o C R d x R d , for some R° > 0. We will prove local existence and uniqueness 
of solutions by a contraction argument in the metric space J- formed by all the 
functions / 6 C([0,T],Vi(R d x R d )) such that the support of f t is contained in 
B R for all t £ [0, T], where R := 2R° and T > is a fixed number to be chosen 
later. Here, we consider the distance in J- given by 

Wi (/,<?):= sup Wxifugt). (28) 
te[o,T] 

Let us define an operator on this space for which a fixed point will be a 
solution to the swarming equation (JTJ) - For / £ J 7 , consider E[f] := Vt/ * p[f]. 
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Then, E[f] satisfies Hypothesis 13.11 (because of the above two Lemmas 13.141 and 
13.151 and the bound (|20[) on Vf) and we can define 



nm) ■■= Tim*/*- 



(29) 



In other words, F[/] is the solution of the swarming equations obtained through 
the method of characteristics, with field E[f] assumed known, and with initial 
condition /o at t = 0. 

Clearly, a fixed point of T is a solution to eq. (QJ on [0, T}. In order for T to 
be well defined, we need to prove that T[f] is again in the space J-, for which 
we need to choose T appropriately. To do this, observe that from eq. (I2"5"j) in 
Lemma 13.141 we have 



\E\ 



Il°°([0,T]x.B b ) - \\^^Wl^(B 2 r 



=:Ci, 



and from point 1 in lemma [3751 



< C 2 



for all P G B R o C R d x R d , and some C?2 > which depends only on a, /3, i?o 
and Ci. Choosing any T < R° jC-x one easily sees that TL^#fo has support 
contained in Br, for all t G [0,T] (recall that we set i? := 2i?°). Then, for each 
t G [0, T], r[/](i) G V\ (R d x R d ), as follows from mass conservation, the support 
of r [/](£) is contained in Br (we just chose T for this to hold), and the function 
t M> T[f](t) is continuous, as shown by Lemma \'3. 121 This shows that the map 
r : T — > T is well defined. 

Let us prove now that this map is contractive (for which we will have to 
restrict again the choice of T) . Take two functions /, g G J 7 , and consider 
r[/],r[g]; we want to show that 



Wi(r[/],r[ ff ]) <cm(f,g) 



(30) 



for some < C < 1 which does not depend on / and g. Using (|2"5|) and (|2U1) . 

Wi(r[/],r[ff])= sup Wxir^f^U^h), (3i) 

and hence we need to estimate the above quantity for each t G [0, T). For 
* G [0, T] , use lemmas |3~TT1 137T1 and [3~T5l to write 



wi(r E[f] #h,r mg] #f ( 



[a] 1 



< 



L=°(supp/ ) 



<C{t) sup \\E[f s ]-E[g s 
se[o,T] 



<C{t)L sup Wi(f s ,g s ) = C(t)LWi(f,g), 
«e[o,T] 
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where C(t) is the function (e c ' 3t — 1)/Cs which appears in eq. (|T7|) . for some 
constant C3 which depends only on a, /3, R, and the Lipschitz constant L of 
VU on B 2R (see eq. pg])). Clearly, 

limCm=0. (32) 

t-!-0 

With (|3T|) . this finally gives 

W 1 (T[f},T[g})<C(T)LW 1 (f,g). 

Taking into account (l3"2"j) , we can additionally choose T small enough so that 
C(T)L < 1. For such T, T is contractive, and this proves that there is a unique 
fixed point of T in T , and hence a unique solution / G J 7 of eq. ([1} . 

Finally, as mass is conserved, by usual arguments one can extend this solu- 
tion as long as the support of the solution remains compact. Since in our case 
the growth of characteristics is bounded (see Lemma T3.2p . one can construct a 
unique global solution satisfying (|2"T1) and (|2"2"|) . □ 



3.4 Stability 

Theorem 3.16. Take a potential U in the conditions of Theorem \3.10\ and fa, 

go measures on R d x R d with compact support, and consider the solutions /, g 
to eq. (HJ) given by Theorem \3.10\ with initial data /o and go, respectively. 

Then, there exists a strictly increasing smooth function r(t) : [0, 00) — > Rq 
with r(0) = 1 depending only on the size of the support of fo and go, such that 



W 1 (f t ,g t )<r(t)Wi(fo,g ), t>0. (33) 

Proof. Fix T > 0, and take R > such that supp/ t and suppgt are contained 
in Br for t s [0,T] (which can be done thanks to theorem 13. 10|) . For t £ [0,T], 
call Lt the Lipschitz constant of T E ^ on Br, and notice that from lemmas 13.81 
and 13.141 we have 

L t <e Clt , te[0,T] (34) 

for some allowed constant C\ > 0. Then we have, using lemmas 13. 1H 13.131 13.71 
and 13131 

W^fugt) = W 1 (T E[f] #f ,T E[g] #9o) 

< ^i(72[/]#/o,T| [s] #/o) + W 1 {T t E[g] #f ,r E[g] #g ) 



< 



'E[f] 



7s [■?] 



L°°(supp/o) 



3 C 2 (i-s) 



Wi(/o,<7o) 



||£[/ s ] -%J|| £ . (Bl0 <fc + Wi(/ ,flb) 

< C 2 Lip 2i ,(VC/) / e c ^WM/ s , 5s ) + e Cl * W^/o 
Jo 



.9o ) 



1G 



Calling C — max{Ci, C2, C2Lip 2i? (V?7)} and multiplying by e ct , 

e^W^fugt) < C f e- Cs W 1 (f s ,g s )ds + W 1 (f ,g ), t e [0,T], 
Jo 

and then by Gronwall's Lemma, 

e- Ct W 1 (f tl g t )<W 1 (f (h go)e C \ te[0,T], 

which proves our result. We point out that the particular rate function r{t) can 
be obtained by carefully looking at the dependencies on time of the constants 
above, leading to double exponentials. □ 

Remark 3.17 (Possible generalizations). As in Remark \3.4\ by assuming more 
restrictive growth properties at infinity of the potential U , we may weaken the 
requirements on the support of the initial data allowing fa with bounded first 
moment for instance. We do not follow this strategy in the present work. 

3.5 Regularity 

If the initial condition for eq. ([1]) is more regular than a general measure on 
M. d x R d one can easily prove that the solution / is also more regular. For 
example, if fo is Lipschitz, then f t is Lipschitz for all t > 0. We will show this 
next. 

Lemma 3.18. Take an integrable function fo : M. d x R d — > [0, +00), with 
compact support, and assume that fo is also Lipschitz. Take also a potential 

UeC 2 (R D ). 

Consider the global solution f to eq. ([1} with initial condition fo given by 
Theorem \3.1(A Then, f t is Lipschitz for all t > 0. 

Proof. Solutions obtained from Theorem 13 . 1 01 have bounded support in velocity 
for all times t > 0, and their fields E = E(t, x) := — VU * p are Lipschitz with 
respect to x. Hence, one can rewrite eq. ((T|) as a general equation of the form 

d t f + div(o/) = 0, 

where a = a(t, x, v) is the expression appearing in the equation, 

a(t, x, v) = (v, E(t, x) + (a — (i |w| 2 )u). 

Then, a is bounded and Lipschitz with respect to x, v on the domain considered 
as the support in velocity is bounded, and classical results show that /( is 
Lipschitz for all t > 0. □ 

4 Well-posedness for General Models 

In this section we want to show that the same results we have obtained in 
the previous section are also valid, with suitable modifications, for much more 
general models than (JlJ . We will start by showing the adaptation of the strategy 
for the Cucker-Smale system and then, we will extend this strategy to more 
general models. 
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4.1 Cucker-Smale Model 



We will prove well-posedness in a slightly more general setting than that of 
the Cucker-Smale model in section [3J being less restrictive on the communi- 
cation rate and the velocity averaging. To be more precise, we shall consider 
£[f](x, v, t) = [(H(x, v)) * /] (x, v, t) as in ©, but for a general H : R d x R d -> 
R d , for which we only assume the following hypotheses: 

Hypothesis 4.1 (Conditions on H). 1. H is locally Lipschitz. 

2. For some C > 0, 

\H(x,v)\<C(l + \x\ + \v\) for all x,v G R d . (35) 

Since the procedure to prove the well-posedness results to © is the same we 
have already applied in the previous section, we will state some of the results 
without proof. First of all, fix T > and let us introduce the system of ODE's 
solved by the characteristics of (J3|): 

V, (36a) 

-&,X,V), (36b) 

where f : [0, T] x R d x R d — >■ R d is any function satisfying the following hypoth- 
esis: 

Hypothesis 4.2 (Conditions on £). 1. £ is continuous on [0, T] x R d x R d , 

5. For some C > 0, 

< C(l + \x\ + \v\), for all t, x, v e [0, T] x R d x R'", and (37) 

5. £ is locally Lipschitz with respect to x and v, i.e., for any compact set 
K C R d x R d i/iere is some > suc/i that 

\Z(t,P 1 )-$(t,P 2 )\<L K \P 1 -P 2 \, *G[0,T], Pi,P 2 eK. (38) 

Under these conditions, we may consider the flow map P| = P|(a;, u) asso- 
ciated to (f36|). defined as the solution to the system (|36]l with initial condition 
(x, v). For ease of notation, we will write the system (f36|) as 

Remark 4.3. Under Hypothesis \4-l\ on H, note that whenever f e C([0, T], Pi(R d x 
R d )) is a given compactly supported measure with supp(/t) C Br* x Brv for all 
te [0, T], the field £[f]=H*f satisfies Hypothesis\f^ 



dt 
d 



-x 



v 
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Definition 4.4 (Notion of Solution). Take H satisfying Hypothesis \4.1\ a mea- 
sure f a £ Vi(R d x R d ), and T £ (0, oo]. We say that a function f : [0,T] -> 
PiQR™ x M. d ) is a solution of the swarming equation^ with initial condition fo 
when: 

1. The field £ = H * / satisfies Hypothesis \4-.2\ 

2. It holds f t = P|#/ . 

Now, an analogue to Lemma 15751 can be stated. We shall state this one and 
the following lemmas for a general £ satisfying Hypothesis (14.21) . 

Lemma 4.5 (Regularity of the characteristic equations). Take T > 0, £ saiis- 
fying Hypothesis (|4.2I) . i? > and t € [0, T]. T/ien there exist constants C and 
L p depending on Lip H (£) and T such that 

|*f(P)|<C for all Pe B R x B R 

and 

|^(Pr) - * 5 (P 2 )| < L p |Px - P 2 | /or a// Pi,P 2 £ B R x B R . 

Lemmas 13.6113.91 are valid as they are presented, taking £ and Hypothesis 
14.21 to play the role of E and Hypothesis 13.11 and making the obvious minor 
modifications on the dependence of the constants. Now we can look at the 
existence of solutions: 

Theorem 4.6 (Existence and uniqueness of measure solutions). Assume H 
satisfies Hypothesis \4-l\ and take fo £ Vi(R d x R d ) compactly supported. Then 
there exists a unique solution f £ C([0,T],Pi(R d x R d )) to equation © in the 
sense of Definition ^ -4] with initial condition fo. Moreover, the solution remains 
compactly supported for all t £ [0,T], i.e., there exist R x and R v depending on 
T, H and the support of fo, such that 



su PP (/ t ) C B R , x B R , for all t £ [0,T]. 



The proof of this result can be done following the same steps as for proving 
Theorem l3.10l Lemmas l3.11l to l5.13l still hold in this situation, and we recombinc 
Lemmas 13 . 141 and 13 . 1 51 in the following result: 



Lemma 4.7. Take H satisfying Hypothesis \4-l\ f £ Vi{ 



with supp(/ t ) C 



Br* x Brv, and £ := £[/] = H * f. Then, for any R > 

LipflXO < Up R+R (H), 
with R := max R X ,R V . Furthermore, if g £ Vi(M. d x R d ) it holds that 



< 



L»(Bfi) 



Up I „. & (H)W 1 (J,g). 



(39) 
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Proof. The first part follows directly from the properties of convolution. For the 
second one, take tt to be an optimal transportation plan between the measures 
/ and g. Then, for any x,v e Br, using that ir has marginals / and g, 

H(x -y,v- u)f(y, u) d(y, u) 



H(x — z, v — w)g(z, w) d(z, w) 

2d 

[H(x — y,v — u) — H(x — z, v — w)] dir(y, u, z, w). 

Taking absolute value, and using that the support of n is contained in the ball 
B k C R id , 

\mx,v)-S\g](x, v)\ 

< / \H(x — y,v — u) — H{x — z,v — w)\ dn(y,u, z,w) 



^ L[ Vr+r( h ) / \(v - z,u-w)\ dir(y,u,z,w) = Lip R+A (H)Wi(f,g). 

□ 

Finally, a stability result also follows using the same steps as in Theorem 

EH 



Theorem 4.8 (Stability in Wi). Assume H satisfies Hypothesis \4-l\ and /q, go £ 
V 1 (R d xR d ) are compactly supported. Consider the solutions /, g to eq. (E|) given 
by Theorem \4-6\ with initial data fo and g , respectively. Then, there exists a 
strictly increasing function r(t) : [0, oo) — 5- Rq with r(0) = 1 depending only 
on H and the size of the support of fo and go, such that 

W 1 {f u9t )<r{t)W 1 {f Q ,go), t>0. (40) 

Remark 4.9 (Evolution of the support in the Cucker-Smale model). In fEJI it 
is shown a sharp bound on the evolution of the support for the kinetic Cucker- 
Smale equation in which H(x,v) = w(x)v. More precisely, it is proved that for 
any given fo 6 Vi(R d x R d ) compactly supported, we have that 

supp(/ t ) C B(x c (0) + mt, R x (t)) x B(m, R v (t)), 

with 

R x (t) < R and R v (t) < R e~ xt 

for some R depending only on Rq = max{_R 2: (0), R v (0)} and A = w(2R). Here, 
m stands for the mean velocity of the system 

v f(t, x, v) dx dv, 
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which is preserved along its evolution. This precise bound on the support and 
the particular choice of H lead to a uniform control in time of the constants 
hip R (H) and L p in the results above, which are now bounded for all times. A 
tedious but straightforward computation leads to a rate r(t) in the stability result 
which is exponentially increasing. Indeed, if we follow the steps of the proof of 
Theorem \3.16\ for the particular case of the Cucker-Smale model we can see 
that, since R x and R v are not increasing with time, the numbers C\ and Ci 
that appear there can be chosen independently of time, whence r(t) shall grow 
at most exponentially. 

Remark 4.10 (Comparison with Literature). As already mentioned above, the 
particular case of the kinetic Cucker-Smale model has already been approached in 
fWjj where the authors give a well-posedness result based on the bounded Lipschitz 
distance. Here, we recover the same result but based on the stability in the 
Wasserstein distance W\, which allows us to obtain sharper constants and rates. 

4.2 General Models 

With the techniques used in the previous sections one can include quite general 
kinetic models in the well-posedness theory. In this section we illustrate this by 
giving a result for a model which includes both the potential interaction and 
self-propulsion effects of section [3J the velocity-averaging effect of section 14.11 
and the more general models above [24l [25] . 

Let us introduce some notation for this section: V c (M. d x R d ) denotes the 
subset of V\ (M d x M. d ) consisting of measures of compact support in W l x M. d , 
and we consider the non- complete metric space A := C([0,T],V c (R d x R d )) 
endowed with the distance W\. On the other hand, we consider the set of 
functions B := C([0, T], Lip, oc (R d x R d ,R d )), which in particular are locally 
Lipschitz with respect to (x,v), uniformly in time. We consider an operator 
%[■] : A — > B and assume the following: 

Hypothesis 4.11 (Hypothesis on a general operator). Take any Rq > and 
f,g £ A such that supp(/t) U supp(gt) C Br for all t £ [0, T], Then for any 
ball Br C K d x R d , there exists a constant C = C(R,Rq) such that 

max \\U[f]-H[g]\\ L ~ (BR) < CW 1 {f,g), 
maxLi PiJ (^[/])<C. 

Associated to this operator, we can consider the following general equation: 

d t f + v-V x f-V v -[U[f]f) = Q. (41) 

Remark 4.12 (Generalization). It is not difficult to see that the choices ~H[f] — 
(a — f3\v\ 2 )v — WU * p and Tl[f] = H * f correspond to (fTJ) and ([3]), respectively, 
and that they satisfy Hypothesis \4-H\ if we assume the hypotheses of Theorems 
\3.1CH and \4-6\ respectively. Moreover, one can cook up an operator of the form: 

H[f] = F A (x, v) + G{x) * p + H(x, v) * / 
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with Fa, G and H given functions satisfying suitable hypotheses, such that the 
kinetic equation (|4ip corresponds to the model (|4]). 

We will additionally require the following: 

Hypothesis 4.13 (Additional constraint on H). Given f € C([0, T], V c (Br )), 
and for any initial condition (X°, V Q ) G M. d x M. d , the following system of ordi- 
nary differential equations has a globally defined solution: 

Tt x = v > ^ 

lv = H[f](t 1 X,V), (42b) 
X(0) = X°, V(0) = V°. (42c) 

Of course, this is a requirement that has to be checked for every particular 
model, and it is difficult to give useful properties of % that imply this and are 
general enough to encompass a range of utile models; therefore, we prefer to 
give a general condition which reduces the problem of existence and stability to 
the simpler one of existence of the characteristics. 

In the above conditions one can follow a completely analogous argument to 
that in the proof of Theorems 13.101 and 13.161 and obtain the following result: 

Theorem 4.14 (Existence, uniqueness and stability of measure solutions for 
a general model). Take an operator %[•] '■ A — > B satisfying Hypotheses \4-H\ 
and \4-l3[ and fo a measure on M rf x M rf with compact support. There exists a 
solution f on [0, +oo) to equation (|41[) with initial condition fo. In addition, 

/GC([0,+4? c (K d xR d )) (43) 

and there is some increasing function R = R(T) such that for all T > 0, 

supp ft C B R(T) C R d x R d for all t e [0, T]. (44) 

This solution is unique among the family of solutions satisfying (1431) and (|44[) . 

Moreover, given any other initial data g 6 P C (M xl ) and g its correspond- 
ing solution, then there exists a strictly increasing function r(t) : [0, oo) — > Rq~ 
with r(0) = 1 depending only on % and the size of the support of fo and go, 
such that 

W 1 (J t ,g t )<r(t)Wi(fo,gol t > 0. 
5 Consequences of Stability 

5.1 iV-Particle approximation and the mean-field limit 

The stability theorems l3.16l and l4.81 or the general version B.Ml give in particular 
a justification of the approximation of this family of models by a finite set of 
particles satisfying a system of ordinary differential equations. We will state 
results for the general model (|4T|). under the conditions on % from section |4~2| 
One can easily check that the following holds: 
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Lemma 5.1 (Particle solutions). Assume H satisfies the conditions of Theorem 
\4-14\ Take N positive numbers mi, . . . , , and consider the following system 
of differential equations: 

x i =v l , i = l,...,N, (45a) 

v i = J2mMf N Kt^i,v i ), i = l,...,N, (45b) 

where f N : [0,T] -)• 7>i(]R d x R d ) is the measure defined by 

N 

ft* ~ ^2 m i S (x i (t),v i (t))- ( 46 ) 
i=l 

If Xi,Vi : [0,T] — > M. d , for i — 1. . . . , N, is a solution to the system (145|) . then 
the function f N is the solution to (|4ip with initial condition 



N 



fo = m i S (xi (0)M0)) ■ ( 4? ) 
i=l 

As a consequence of the stability in W\ , we have an alternative method to 
derive the kinetic equations (fTJ) , Q or (T4H) , based on the convergence of particle 
approximations, other than the formal BBGKY hierarchy in [31 [7]. 

Corollary 5.2 (Convergence of the particle method). Given fo G Vi(M. d x M. d ) 
compactly supported and f-L satisfying the conditions of Theorem \4-14\ take a 
sequence of f^ of measures of the form (|47[) (with mi, 2^(0) and Uj(0) possibly 
varying with N), in such a way that 

lim Wi(/ ",/o) = 0. 

N — >QC 

Consider f^ given by (1461) . where Xi(t) andvi(t) are the solution to system (|45p 
with initial conditions Xi(0), Vi(0). Then, 

lim W 1 (f t N J t ) = 0, 

N— >oo 

for all t > 0, where f — f(t, x, v) is the unique measure solution to eq. (|4ip 
with initial data fo. 

5.2 Hydrodynamic limit 

We state our hydrodynamic limit result for eq. ([1} . If we look for solutions of 
([1} of the form 

f(t, x, v) = p(t, x) 6(v - u(t, x)) (48) 

for some functions p, u : [0, T] x R d — > M, one formally obtains that p and u 
should satisfy the following equations: 

d t p + div x {pu) =0, (49a) 
d t u + (u ■ V)u = u(a-/3 \u\ 2 ) -VU*p. (49b) 
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This is made precise by the following result whose existence part was already 
obtained in [7]: 

Lemma 5.3 (Uniqueness for Hydrodynamic Solutions). Take a potential U G 
C 2 (E d ) and assume that there exists a smooth solution (p, u) with initial data 
(po,Uo) to the system (|49p defined on the interval [0, T]. Then, if we define 
f : [0,+oo) -> V\ {M. d x M. d ) by 

/ f(t,x,v)<f)(x,v)dxdv = / <j){x, u(t, x)) pit,, x) dx (50) 

JR d xR d JR d 

for any test function cj) G C^,(M. d x R d ) 7 then f is the unique solution to ((TJ) 
obtained from Theorem \3.10\ with initial condition fo = PqS(v — uq). 

As a direct consequence of Lemma 15.31 and the stability result in Theorem 
13.161 we get the following result. 

Corollary 5.4 (Local-in-time Stability of Hydrodynamics). Take a potential 
U G C 2 (M. d ) and assume that there exists a smooth solution (p,u) with initial 
data (pq,uq) to the system (|49l) defined on the interval [0,T]. Let us consider a 
sequence of initial data /q G ViQSL d X M. d ) such that 

lim Wxif^poSiv-Uo)) =0. 

k^oc 

Consider the solution f k to the swarming eq. ([TJ) with initial data f§. Then, 

lim W 1 (tfj t )=0, 

for all t £ [0, T] with f(t, x, v) — p(t, x) S(v — u(t, x)). 

Acknowledgments 

The authors acknowledge support from the project MTM2008-06349-C03-03 
DGI-MCI (Spain) and 2009-SGR-345 from AGAUR-Generalitat de Catalunya. 
This work was developed at the CRM-Barcelona during the thematic program 
in "Mathematical Biology" in 2009. 

References 

[1] A. Barbaro, K. Taylor, P. F. Trethewey, L. Youseff and B. Birnir, Discrete 
and continuous models of the dynamics of pelagic fish: application to the 
capelin, preprint. 

[2] B. Birnir, An ODE model of the motion of pelagic fish, J. Stat. Phys. 128 
(2007) 535-568. 

[3] N. N. Bogoliubov, Problems of a dynamical theory in statistical physics, 
Studies in Statistical Mechanics, Vol. I, (Interscience, New York, 1962). 



24 



W. Braun and K. Hepp, The Vlasov Dynamics and Its Fluctuations in 
the 1/N Limit of Interacting Classical Particles, Commun. Math. Phys. 
56 (1977) 101-113. 

E. Bonabeau, M. Dorigo and G. Theraulaz, Swarm Intelligence: From 
Natural to Artificial Systems, (Oxford University Press, New York, 1999). 

S. Camazine, J.-L. Deneubourg, N.R. Franks, J. Sneyd, G. Theraulaz and 

E. Bonabeau, Self- Organization in Biological Systems, (Princeton Univer- 
sity Press, 2003). 

J. A. Carrillo, M. R. D'Orsogna and V. Panferov, Double milling in self- 
propelled swarms from kinetic theory, Kinetic and Related Models 2 (2009) 
363-378. 

J. A. Carrillo, M. Fornasier, J. Rosado and G. Toscani, Asymptotic Flock- 
ing Dynamics for the kinetic Cucker-Smale model, preprint UAB, . 

J. A. Carrillo and G. Toscani, Contractive probability metrics and asymp- 
totic behavior of dissipative kinetic equations, Riv. Mat. Univ. Parma 6 
(2007) 75-198. 

Y.-L. Chuang, M. R. D'Orsogna, D. Marthaler, A. L. Bertozzi and L. 
Chayes, State transitions and the continuum limit for a 2D interacting, 
self-propelled particle system, Physica D 232 (2007) 33-47. 

I.D. Couzin, J. Krause, N.R. Franks and S.A. Levin, Effective leadership 
and decision making in animal groups on the move, Nature 433 (2005) 
513-516. 

F. Cucker and S. Smale, On the mathematics of emergence, Japan. J. 
Math. 2 (2007) 197-227. 

F. Cucker and S. Smale, Emergent behavior in flocks, IEEE Trans. Au- 
tomat. Control 52 (2007) 852-862. 

R. Dobrushin, Vlasov equations, Fund. Anal. Appl. 13 (1979) 115-123. 

M.R. D'Orsogna, Y.-L. Chuang, A. L. Bertozzi, L. Chayes, Self-propelled 
particles with soft-core interactions: patterns, stability, and collapse, 
Phys. Rev. Lett. 96 (2006) 104302-1/4. 

R. Eftimie, G. de Vries, M.A. Lewis, Complex spatial group patterns 
result from different animal communication mechanisms, Proceedings of 
the National Academy of Sciences 104 (2007) 6974-6979. 

F. Golse, The Mean- Field Limit for the Dynamics of Large Particle Sys- 
tems, Journees equations aux derivees partielles 9 (2003) 1-47. 

G. Gregoire and H. Chate, Onset of collective and cohesive motion, Phy. 
Rev. Lett. 92 (2004) 025702-1/4. 



25 



[19] S.-Y. Ha and J.-G. Liu, A simple proof of the Cucker-Smalc flocking dy- 
namics and mean- field limit, to appear in Comm. Math. Sci. 

[20] S.-Y. Ha and E. Tadmor, From particle to kinetic and hydrodynamic 
descriptions of flocking, Kinetic and Related Models 1 (2008) 415-435. 

[21] M. Hauray and P.-E. Jabin, iV-particles approximation of the Vlasov equa- 
tions with singular potential, Arch. Ration. Mech. Anal. 183 (2007) 489- 
524. 

[22] A.L. Koch, D. White, The social lifestyle of myxobacteria, Bioessays 20 
(1998) 1030-1038. 

[23] H. Levine, W.-J. Rappel and I. Cohen, Self-organization in systems of 
self-propelled particles, Phys. Rev. E 63 (2000) 017101-1/4. 

[24] Y.X. Li, R. Lukeman and L. Edclstein-Keshet, Minimal mechanisms for 
school formation in self-propelled particles, Physica D 237 (2008) 699- 
720. 

[25] R. Lukeman, Y.X. Li and L. Edelstein-Keshet, A conceptual model for 
milling formations in biological aggregates, Bull Math Biol. 71 (2008) 
352-382. 

[26] A. Mogilner, L. Edelstein-Keshet, L. Bent and A. Spiros, Mutual interac- 
tions, potentials, and individual distance in a social aggregation, J. Math. 
Biol. 47 (2003) 353-389. 

[27] H. Neunzert, The Vlasov equation as a limit of Hamiltonian classical me- 
chanical systems of interacting particles, Trans. Fluid Dynamics 18 (1977) 
663-678. 

[28] H. Neunzert, An introduction to the nonlinear Boltzmann- Vlasov equa- 
tion, in Kinetic theories and the Boltzmann equation (Montecatini, 1981), 
pp. 60-110, Lecture Notes in Math. 1048, (Springer, Berlin, 1984). 

[29] J. Parrish, L. Edelstein-Keshet, Complexity, pattern, and evolutionary 
trade-offs in animal aggregation, Science 294 (1999) 99-101. 

[30] L. Perea, G. Gomez, and P. Elosegui, Extension of the Cuckcr-Smale con- 
trol law to space flight formations, AIAA Journal of Guidance, Control, 
and Dynamics 32 (2009) 527-537. 

[31] H. Spohn, Kinetic equations from Hamiltonian dynamics: Markovian lim- 
its, Rev. Modern Phys. 52 (1980) 569-615. 

[32] H. Spohn, Large scale dynamics of interacting particles, Texts and Mono- 
graphs in Physics, Springer 1991. 



26 



[33] CM. Topaz and A.L. Bcrtozzi, Swarming patterns in a two-dimensional 
kinematic model for biological groups, SIAM J. Appl. Math. 65 (2004) 
152-174. 

[34] CM. Topaz, A.L. Bertozzi and M.A. Lewis, A nonlocal continuum model 
for biological aggregation, Bulletin of Mathematical Biology 68 (2006) 
1601-1623. 

[35] T. Vicsek, A. Czirok, E. Ben- Jacob, I. Cohen and O. Shochet, Novel type 
of phase transition in a system of self-driven particles, Phys. Rev. Lett. 75 
(1995) 1226-1229. 

[36] C. Villani, Topics in optimal transportation, Graduate Studies in Mathe- 
matics Vol. 58, (Amer. Math. Soc, Providence, 2003). 



27 



